Anomalous spin density distribution on oxygen and Ru in Cai 5Sro.5R.uO4: 
A polarised neutron diffraction study 
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By means of polarized neutron diffraction in a magnetic field of 7.0 T at f .6 K an anomalously large 
magnetization density is observed on the in-plane oxygen in Cai.sSro.sRuCU. Field-induced moments 
of different ions are determined by refinement on the flipping ratios, yielding /j,r u = 0.346(11) /j,b, 
Hoi = 0.076(6) fj,B and po2 = 0.009(6) fj,s- The moment on the oxygen arises from the strong 
hybridization between the Ru-4d and 0-2p orbitals. The maximum entropy magnetization density 
reconstruction reveals a strongly anisotropic density at the Ru site, consistent with the distribution 
of the xy (t2g band) rf-orbitals. 

PACS numbers: 75.25.+Z, 75.50.-y, 75.50.Cc, 75.50.Gg 



Layered perovskite ruthenates have attracted consid- 
erable interest because of their magnetic properties and 
in particular because of the recent discovery of supercon- 
ductivity in Sr2RuC>4 It is now widely accepted that 
superconductivity in this compound is unconventional 
and the pairing of triplet p-wave symmetry. Sr2Ru04 
does not exhibit any magnetic ordering, but among the 
related compounds with a Ru 4+ -ion various magnetic 
structures are observed : ferromagnetic order occurs in 
the metallic perovskite SrRuC>3 Q, antiferromagnetic 
one in insulating Ca2RuC>4 ||] and incommensurate order 
in Sr2Ruo.91Tio.09O4 {§. 

According to the band-structure calculations, the 
Fermi-surface in Sr2Ru04 is formed by three bands as- 
sociated with the three t2 g Ru orbitals, xy,yz and xz ||. 
Since 4d i° ns generally have more extended d orbitals 
than the corresponding 3d ions, the ^rf-oxides exhibit 
larger overlap and hybridization between the transition 
metal and O 2p orbitals. As a result magnetism becomes 
more itinerant in character and there is stronger interplay 
between structural degrees of freedom, and the magnetic 
and electronic properties. Another important point is 
that in 3d transition-metal systems, the strong crystal 
field interaction dominates over the spin-orbit interac- 
tion and consequently the orbital moment is quenched, 
fiL = 0, whereas, in layered ruthenates, much stronger 
spin-orbit effects may be anticipated, which should man- 
ifest themselves in the radial part of the Ru form factor. 

The substitution of Sr by Ca in Sr 2 Ru04 leads to a 
complex phase diagram || |7| . Already small concentra- 
tions of Ca give rise to appearance of a structural dis- 
tortion characterized by the rotation of RuOg-octahedra 
around the c-axis, which is accompanied by a strong in- 
crease of the low temperature magnetic susceptibility. 
For Sr-concentrations higher than 0.5 in C^-rrSr^RuG^ 
additional distortions are observed characterized by a tilt 
of the octahedra around an in-plane axis. In contrast to 



the rotational distortion, the increase of the tilt leads to 
a decrease of the low temperature magnetic susceptibil- 
ity. The critical point x=0.5, where the tilt distortion 
appears, corresponds to the maximum of the low tem- 
perature susceptibility. The susceptibility at this point 
is about 100 times larger than that of pure Sr 2 Ru04 [||. 
We choose this composition for our study, since the large 
susceptibility is favorable for the polarized neutron mea- 
surements. 

Polarized neutron diffraction (PND) is a uniquely pow- 
erful tool to study the magnetization densities in crystals, 
due to the fact that it provides direct information about 
the 3D distribution of the magnetization throughout the 
unit cell. In favorable cases, PND further allows to de- 
termine the symmetry of occupied orbitals. The aim of 
the present work is two-fold : i) to check the presence 
of unpaired electron density on the ligands 01 (and 02) 
which the theory predicts to be anomalously high, and 
ii) to determine the radial part of the Ru magnetic form 
factor and possibly to establish the symmetry of Ru 4d 
orbitals occupied by unpaired electrons. 

Polarized neutron measurements were performed on 
the two-axis lifting-counter diffractometer 5C1 at the 
ORPHEE reactor, LLB CEA/Saclay, using a wavelength 
A =0.845 A (Heusler alloy monochromator, polarization 
of the incident neutron beam Po = 0.91). Higher-order 
contamination was suppressed to less than 0.01% by 
means of erbium filters. The flipping ratio (ratio between 
neutron spin-up and spin-down intensity) measurements 
provided the experimental data used in the subsequent 
magnetic moment refinement and in the magnetization 
reconstruction. The programs MAGLSQ and SORGAM 
of the Cambridge Crystallography Subroutine Library |J 
were used for least squares refinement on the flipping ra- 
tios and for the calculation of magnetic amplitudes, re- 
spectively. 

A single crystal of Cai.5Sro.5Ru04 was obtained by 
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FIG. 1: Magnetic susceptibility of Ca1.5Sro.5O4 obtained from 
the polarized neutron intensities of the (004) reflection with 
the field parallel to the [110] direction. 



a floating zone technique g]. According to Ref. 
Cai.5Sro.5Ru04 crystallizes in the tetragonal IAi/acd 
structure, for which the oxygen octahedra are rotated 
around the c axis with an opposite phase of rotation for 
the planes separated by 12.5 A. The structural parame- 
ters of the Cai.5Sro.5Ru04 crystal used in our experiment 
can be summarized as follows: a=6=5.3A, c=25.0 A, Ru 
8a position, Ca(Sr) 16d z=0.549, 01 16/ x=0.193, 02 
16d z=0.457 (see § for details). 

The temperature dependence of the magnetic suscep- 
tibility x °f our Cai.5Sro.5Ru04 crystal was measured 
directly by means of PND in magnetic external fields of 
1 T, 4 T and 7 T. For this purpose, the spin-up and 
spin-down intensities of the (004) reflection were 

collected as a function of temperature, and the difference 
I + — I~ plotted versus temperature (see Fig{l]). This dif- 
ference, I + — I~ = 2Po\HFn, where Fn is the nuclear 
structure factor for the same reflection and Pq the po- 
larization of the incoming beam, is proportional to the 
Fourier component (004) of the wave-vector dependent 
susceptibility of the crystal. As seen from Fig.l and 
magnetization data, the susceptibility of Cai.sSro.5Ru04 
increases strongly at low temperatures but the crystal 
remains paramagnetic down to 1.6 K; at 7 T the magne- 
tization is near saturation. 

Overall, two sets of polarized neutron flipping ratios 
were measured at 1.6 K. Regarding the first, a verti- 
cal magnetic field of 7 Tesla was applied parallel to the 
[110]-axis of the crystal and a total of 94 flipping ratios 
with sin#/A < 0.6 A -1 . This yielded 53 independent 
reflections of which only a subset of 13 reflections has a 
pure oxygen contribution. (Note that oxygen atoms con- 
tribute to all the reflections since they occupy low sym- 
metry 1 6f and 16e positions.) The second set consisted 
of 97 (34 independent and among them 13 pure oxygen 
reflections) flipping ratios also measured at 1.6 K, but 
with the magnetic field along the [001] direction. 

The magnetic amplitudes of all measured reflections 
were obtained using the SORGAM |j| program, except 
for those containing a pure oxygen contribution. They 
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FIG. 2: Magnetic amplitudes of Ca1.5Sro.5O4 at 1.6 K and 
7 T. Circles represent mixed reflections, normalized by Ru 
geometrical structure factor. Squares in the inset are the 
pure oxygen reflections normalized by the geometrical struc- 
ture factor of 01. The solid line is the magnetic amplitude 
calculated for the Ru form factor in the dipole approxima- 
tion. Thin lines show the decomposition of the latter into its 
components proportional to fiL and (fis + ^l) respectively. 



are shown in Fig. || together with the theoretical mag- 
netic form factor of Ru. (For convenience, the ampli- 
tudes were normalized by the geometrical structure fac- 
tor of the Ru.) As seen from Fig. 2, the majority of 
measured magnetic amplitudes lie reasonably close to the 
theoretical curve. There are, however, a large number of 
reflections, especially those at low Q, for which the dis- 
agreement between experimental and theoretical values is 
very significant. A simple analysis shows that these are 
the reflections with a strong contribution of the 01 sub- 
lattice. Those amplitudes of 16 reflections containing a 
pure oxygen contribution (h+k=2n+l) were normalized 
by the geometrical structure factor of 01. The ampli- 
tudes calculated in such a way should be proportional to 
the oxygen form factor, if one neglects the aspherical part 
of Ru magnetic density. The comparison of experimental 
results (see inset in Fig.||) with the solid line represent- 
ing the theoretical form factor of oxygen shows that an 
additional moment of the order of 0.075 hb is present at 
the 01 site. 

Further analysis of the measured flipping ratios was 
performed by refinement. The refinement of the first data 
set using the tetragonal group IAx/acd demonstrated 
clearly that there is a very large moment induced both 
on Ru and on 01. The moments induced by the field 
of 7 T were obtained as = 0.350(15) /is, fioi — 
0.076(8) /i S and /i 02 = 0.014(8) fi B , with X 2 = 1-87 for 
94 reflections. For comparison a model in which the oxy- 
gen moment is neglected gives [1r u — 0.33(15) [Ib and 
X 2 = 10.2. The refinement performed using only 13 pure 
oxygen reflections yields fioi = 0.072(12). The refine- 
ment of the second data set with the applied magnetic 
field parallel to the [001] direction gives a very similar 
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result [Xr u = 0.342(15) /is, /ioi = 0.076(8) [Ib and \ioi 
= 0.009(6) hb with x 2 = 1-98 for 97 reflections. 

These refinements show unambiguously that apart 
from the relatively large magnetic density induced at the 
Ru position a huge amount of density exists at the 01 
site. This moment amounts to more than 20 % of that of 
Ru. Taking the multiplicity of the 01 site into account, 
one can conclude that nearly one third of the total mag- 
netic density in the cell is transferred from Ru to 01. 
By contrast, the refinement does not show any evidence 
for the presence of any magnetic moment neither at the 
apical oxygen 02 nor at the Sr sites. 

In the least-squares analysis, the Ru magnetic form 
factor in the dipole-approximation was used 

mufRu{Q) = (us + fJ>L){jo(Q)) + fJ>L(j2(Q)). Here, 
US and are the spin and the orbital components of the 
total magnetic moment hr u and (jo(Q)) and (72 (Q)) are 
radial functions. Upon increasing Q, (jo(Q)} decreases 
whereas (72 (Q)) has a maximum at Q=0.3 A" 1 . We 
used in the refinement the radial functions of the Ru 1+ 
state, since those of higher ionic states of Ru are not 
available. At first sight, the utilization of the radial 
integral of the ionic state for ruthenates, being itinerant 
in nature and often mixed strongly with conduction 
electrons, is not well justified. We remind the reader, 
however, that a similar problem exists in actinide 5f 
compounds also characterized by a strong hybridization. 
For actinides, it has been found that the radial part of 
the wave functions does not change appreciably through 
the actinide series and for different ionic states. On 
the contrary the mixing coefficient C2 = ^l/(^s + Hl) 
depends strongly on the hybridization ]ll[ | . The validity 
of this approach for J^d compounds remains to be verified 
but the description of the Ru form factor by mixing of 
O'o(0)) and (72 (Q)) of the Ru 1+ state can be considered 
as a good approximation. The value of coefficient C2 
will strongly depend on the ionic state and hybridization 
of Ru, and as in actinides, it can be obtained from 
refinement against the flipping ratios. In our instance, 
the orbital component of the Ru moment was included 
and was found to be fiL = 0.140(15) fig- Neglect of 
the orbital part of the Ru form factor deteriorated 
the quality of the fit considerably, from \ 2 = 1.98 to 
X 2 = 2.88, thereby indicating that the orbital part of 
the moment exists and should be taken into account. 
The decomposition of the Ru form factor into the radial 
functions (jo(Q)) and (72 (Q)) is shown in Fig.|[ 

The results of our refinements show a significant spin 
density at the 01 site and it is useful to check to what 
extent this can be substantiated by a model-free analysis 
of our data, obtained for instance by reconstructing the 
3D magnetization distribution using the maximum en- 
tropy method (MEM). This method has been shown to 
give much more reliable results than conventional Fourier 
syntheses, by considerably reducing both noise and trun- 
cation effects |lCj . In order to carry out the MEM recon- 
structions, the magnetic structure factors of both data 
sets, related to two different orientations of the sample, 



were merged into 70 independent flipping ratios. The 
magnetization density distribution was discretized into 
48x48x240 sections along a,b and c respectively. Then 
it was reconstructed using a conventional uniform (flat) 
density prior p0[ . Such a procedure is biased against the 
creation of any magnetic density in the unit cell. In spite 
of the negative bias of the MEM procedure the ensuing 
reconstruction revealed the presence of strong peaks at 
the positions of Ru and weaker but well pronounced ones 
at the 01 positions. Numerical integrations of the mag- 
netization density over the Ru and the 01 positions yield 
fJ>Ru = 0.307 jUb, /xoi = 0. 030^_b. Because of the nega- 
tive bias mentioned above, the MEM values are slightly 
smaller than those from the refinement. Moreover, sec- 
tions through the Ru site revealed a strongly anisotropic 
magnetization density, while sections through the Sr and 
the 02 sites did not detect any magnetic density above 
the 3 x 1Q~ 3 hb level. To check the robustnest of the den- 
sity reconstruction against possible spurious features, the 
MEM procedure was repeated using a non-uniform (ad- 
verse) density prior |l2L in which all the magnetisa- 
tion in the crystal is spherically concentrated around the 
Ru sites. This non-uniform prior was calculated analyt- 
ically using the spherical magnetic density of Ru ions at 
the 8a position via the moment and radial integrals val- 
ues obtained from the least-squares refinement mentioned 
previously. Such a procedure considerably enhances the 
bias against both the creation of magnetic density at ei- 
ther oxygen sites and the asphericity of Ru density, as 
compared to the uniform prior strategy. Nevertheless, 
both features of interest, the oxygen density and the as- 
phericity of the Ru distribution, survive the acid test of 
an unfavorable non-uniform prior. The final results us- 
ing the non-uniform prior on the basis of all measured 
flipping ratios are shown in Fig||. 

Based on the band-structure calculations a strong hy- 
bridization between Ru-4d and the 0-2p orbitals was 
postulated for the ruthenates in general ||. In conse- 
quence a significant oxygen contribution was introduced 
into the Stoner interaction parameter for Sr2Ru04. Since 
the oxygen may not be polarized in an antiferromagnetic 
fluctuation, the Stoner interaction is then q-dependent. 
Taking this Stoner-interaction into account the magnetic 
excitation spectrum for Sr2Ru04 has been calculated and 
used as an input for many theories to explain supercon- 
ductivity in this compound ^ |TJ, |l5|, [l6| . Our finding of 
a large spin-density on the oxygen-site validates the basis 
of these procedures. So far, there are no calculations of 
the spin-density induced by an external field available, 
but the magnetization in two rather distinct ferromag- 
netically ordered ruthenates was calculated ||, SrRuO,3 
and Sr2RuYOg. In both compounds one third of the to- 
tal spin desity was found on the oxygen sites, in perfect 
agreement with the value we observe for the paramag- 
netic compound. A comparable amount of transferred 
magnetization has so far not been observed in any oxide, 
which underlines this outstanding property of the ruthen- 
ates compared to the better studied 3d-compounds. For 
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FIG. 3: Maximum entropy reconstruction of the magnetization distribution in Cai.sSro.sRuCU at 1.6 K and 7 T. The maps 
A and B show the sections passing through the Ru02-plane at z=| and map C shows the section perpendicular to the plane 
through the Ru-position. The edge lengths of the maps correspond to the tetragonal lattice constant of a=5.34A for maps A 
and B and to a and a quarter of c for map C respectively. Note that in space-group I4i/acd a Ru is situated at and an 

Ol-site lies at (0. 19,0. 44, |). The section in A) relates to the adverse 3D non-uniform prior density used to compute B) and C), 
which biases the MEM reconstruction against any magnetization at the O sites as well as against an aspherical magnetization 
at the Ru site. The contour step is 0.03/xs A~ 3 . 



example, PND found an induced moment of only a few 
percent in an Fe-garnet []l7f as well as in two manganates 

m 

The anisotropic density distribution at the Ru is ex- 
tended along the diagonals of the original planar per- 
ovskite lattice (i.e. in the direction 45 degrees to the 
bonds) and it is flattened along the c-direction, see Fig. 
3. This anisotropy suggests that the main part of the 
spin-density arrises from the two-dimensional 7-band re- 
lated to the planar xy-orbitals. In general the 7-band is 
considered as being the driving element for a ferromag- 
netic instability in the 214-ruthenates, since there is a 
van Hove singularity slightly above the Fermi-level [^). 
Fang and Terakura have analyzed the phase diagram of 
Ca2- K Sr 2: Ru04 by ab-initio calculations based on the lo- 
cal density approximation jnj]. They observe that the 
rotation of the octahedra flattens and slightly lowers the 
7-band which enhances the influence of the van Hove 



singularity as well as the instability towards ferromag- 
netism. The anisotropy observed for Cai.5Sro.sRu04, a 
sample very close to ferromagnetic order, is fully consis- 
tent with this explanation. However, it does not agree 
with the result of the LDA plus U calculation by Anisi- 
mov et al. for that particular concentration po| . 

In conclusion, we have unambiguously evidenced a sig- 
nificant induced spin density at the positions of Ru hr u 
= 0.350(15) and anomalously high spin density at the in- 
plane oxygen 01 \io\ = 0.076(8) [Ib in Cai.sSro.sRuO^ 
In contrast, neither the least-squares refinement nor the 
model-free MEM procedure show a significant magnetic 
moment at the sites of Sr and apical oxygen 02. The 
large magnitude of the moment transferred to 01 is con- 
sistent with theorical predictions ||. 
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ence and Technology Corporation, and by Deutsche 
Forschungsgemeinschaft through SFB 608. 
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